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Abstract 

A new event generator is proposed for two processes e + e _ — ► vv^ 
and e + e~ — ► Z/P77 where vv includes all the neutrino species. The 
exact matrix elements of single- and double-photon emission, gen- 
erated by the GRACE system, are convoluted with the QED parton 
shower (QEDPS) to deal with the initial state radiations (ISR). It is 
pointed out that a careful treatment is required to avoid the dou- 
ble counting of the radiative photons between the matrix elements 
and the ISR part. A detailed comparison of with the 0(a) 

calculations and other similar Monte Carlo generators is discussed on 
the total cross section and on various distributions. It is also examined 
how the possible effects of the anomalous triple-gauge-boson couplings 
can be observed. 
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1 Introduction 



The processes of single- and double-photon emission associated with a large 
missing energy are very important for the precision test of the standard model 
and the search for any new (invisible) particles beyond the standard model. 
In fact the number of the light-neutrino generations has been measured us- 
ing this process on the Z-pole at the LEP experiments]!]]. The efforts to 
find new invisible particles, such as a neutral (lightest) SUSY particle or a 
sequential heavy neutrino, are now continued by the LEP-2 experiments at 
higher energies [@]. In this case, however, the neutrino pair-productions with 
hard-photon emission are the main background for the investigation. With 
the increase of the accumulated luminosity, the cross sections of these pro- 
cesses must require more precision with an uncertainty of (9(1%). In addition 
since e + e~ — > VeV^i^) includes the triple (quartic) gauge couplings, there is 
a possibility to find some new effects beyond the standard model on these 
rather poor-tested couplings. 

There are some methods to estimate the cross sections including the 
higher order QED corrections. The simplest one would be to combine the 
matrix elements of e + e~ vv with some tool for the initial-state radia- 
tions(ISR). The hard photon(s) may be supplied, for instance, by a parton 
shower, which can treat the photons with a finite transverse moment um(p'r)- 
When this method is applied to the center-of-mass system(CMS) energies 
above the Z-boson mass, however, the leading logarithmic (LL) approxima- 
tion is not sufficient. This is because the colliding energy after the photon 
emission can match the Z-pole('Radiative Return'), and the probability to 
emit high-energy photon(s) is largely enhanced. When this happens the LL 
approximation is not successful in describing the spectrum around the Z 
pole. To improve the precision of the calculations the non-LL terms for ISR 
must be taken into account as well as the LL terms. However, there exist no 
complete calculations up to the non-LL order. 

Another possible way is to make use of the matrix elements of e + e~ — ► 
^i?7(7)Q together with the ISR tool. Diagrams of the single(double)-photon 
emission are shown in Fig.l (Fig. 2), respectively. The visible and invisible 
hard-photon(s) are fed by the matrix elements as well as by the ISR tools. 

1 The complete matrix elements of the e + e~ — > v^v^ process and their compact 
approximated form were first given in Ref.jjj). 
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Here it is crucial to avoid the double counting in the radiations coming from 
these two different sources. The conventional QED structure function is not 
suitable for this purpose because it is obtained after the integration over all 
allowed phase-space of the emitted photons except for their energy. Thus the 
double counting necessarily happens and cannot be excluded. 

There exists no full electroweak one-loop calculation for e + e~ — > vv^. 
Only the limited corrections, one-loop QED corrections and self-energy cor- 
rections to the Z-exchange diagrams [|, [l| are available. Though there 
is no one-loop calculation at all for the v e v e ^j process, the ISR tools might 
cover the most important corrections and allows us to give a realistic event 
generator which can meet the required theoretical uncertainty of (9(1%). 

The KORALZ Monte Carlo program || exploits the exact matrix elements 
of the single-photon processes with and v T production(an approxima- 
tion is made for u e ) together with the exclusive exponentiation of YFSJ7] 
as the ISR tool. Another event generator is NUNUGVP||, which uses the exact 
matrix elements of the single- and multi(up to three)-photon emission(for 
all kinds of neutrinos) generated by the ALPHA algorithmic and the pr- 
depending structure function IfLOU for ISR. Both programs include anomalous 
triple-gauge-boson coupling(TGC) in the package. 

In this report a new event generator grcz/z/7 is proposed. It adopts the 
exact matrix elements for e + e~ — > 1^7(7) generated by means of the GRACE 



system |TT[] combined with QEDPS|T2| for ISR. The advantages of these pack- 
ages are: 

• The exact matrix elements up to the double-photon emission, includ- 
ing the v e process, are used. Double-photon emission is practically 
sufficient for experimental analysis. 

• QEDPS keeps the complete kinematics for the emitted photons and vir- 
tual electrons before collisions. It allows a more flexible treatment of 
the ISR effects in avoiding the double-counting. 

In section 2, the calculation method, particularly how to connect QEDPS to 
the radiative processes without double counting, is explained. The numerical 
results of single- and double-photon emission are presented in section 3. A 
detailed comparison among grcz/i/7, the exact 0(a) calculations, KORALZ 
and NUNUGPV is also given in this section. The effects of the anomalous TGC 
in the single-photon events are discussed in section 4. The conclusions are 
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summarized in section 5. Appendix describes how to obtain and use the 
grci/wy program. 



2 Calculation Method 

The exact matrix-elements of all the processes are produced by an automatic 
calculation system GRACE. It generates the FORTRAN routines needed to 
calculate the amplitudes numerically based on the helicity-amplitude for- 
malism. 

For ISR the parton shower algorithm for QED, QEDPS, is used. In this 
algorithm the Altarelli-Parisi equation is solved in the LL approximation |T4] 
using the Monte Carlo method. The details of this method can be found in 
Ref . |T2[ . Here we recall that the algorithm can maintain the exact kinematics 



during the evolution of an electron. Let us consider the branching process 

e~(x, K 2 ) -> e-(xy, K 2 ) + 7 (z(l - y), Q 2 Q ), 
where the parent electron, whose virtuality is K 2 , has the momentum p - 



(E, QtiPz), E = ^p 2 z — K 2 and p z = xp* with p* being the momentum of the 
system in the infinite-momentum-frame. A cutoff mass Qq is introduced for 
the photon. Then the momenta of the daughter particles, p\ for e~ and P2 
for 7, can be expressed as 

Pi = (Ei,1s. T ,yp s ), 

P2 = (E 2 ,-k T ,(l-y)p z ), 

where 



Ei = Jy 2 p 2 z + k 2 T -K 2 . 



E 2 = ^{l-yfpl + ^ + Ql 

Assuming p* — > oo, we have 

- K 2 = -K 2 /y + Q 2 /(l -y) + k 2 T /(y(l - y)), 

which determines from y, K 2 and K\. Hence one can obtain the distribu- 
tion of the transverse momentum of the emitted photons, k^. Moreover this 
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enables us to find the virtuality of the electron, —K 2 , as well as its energy, 
just before the collision. 

The phase-space integration of the matrix element squared is carried out 
numerically by BASES (T^j using an adaptive Monte Carlo method. The inte- 
gration is done in the five(eight)-dimensional kinematical phase-space of the 
three (four)-body final state. The evolution of the initial-state electrons by 
QEDPS is treated independently of BASES. 

Since the QEDPS can provide a complete kinematical information about 
the emitted photons and the virtual electrons, it is easy to distinguish the 
photons from the matrix-elements and those from the parton shower. Let us 
divide the full one-photon phase space f2 7 into two regions: One is for the 
visible region, Q v , and the other is for the invisible region, Q±. These are 
defined as 



where E 1 is the energy of the photon and 9 1 is its polar angle in degree 
from the beam axis. The minimum energy and angle, Eq and 6q, must be 
given by the experimental condition. First the photons from QEDPS, 7 ps 's, 
are restricted to go inside of fl± and the photon from the matrix-elements, 
7 ME , is emitted into Q v . Since there is no overlap between these two regions 
no double-counting occurs. In addition the ordering of the electron virtu- 
ality is also required. During the evolution of an electron the virtuality is 
monotonically increasing, which is realized naturally in the QEDPS algorithm. 
A further condition must be imposed on the virtuality of the electron in the 
matrix-elements after emitting the photon: It should be greater than the 
virtuality of the electron in the last stage of QEDPS. 

The configuration mentioned here cannot cover all the configurations for 
the photon emission. When n photons are emitted from one electron line, 
the first (n — 1) photons are fed by QEDPS, including the statistical factor of 
l/(n — 1)!, and only the last photon is supplied by the matrix element. This 
configuration represents the case where the first (n — 1) 7 ps 's are invisible 
but only the last one, 7me, is visible. Another case also exists. One of the 
first (n — 1) 7 ps 's goes into Q v while the other (n — 2) 7 ps 's and 7 M e escape 
into Q±. Though 7me is invisible in this case, the results are free from any 
infrared divergence thanks to the virtuality cut Q introduced above. The 




8ry) \E~ > Eq, 9q < 9j < 180 — 6q}, 



(1) 
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minimum energy of 7me is determined from the minimum virtuality imposed 
by the QEDPS evolution. If there is no 'Radiative Return' the contribution 
of this second configuration is negligible small. However, for the process 
e + e~ — > 1/1/7(7) above the Z° threshold its contribution to the total cross 
section becomes sizeable, typically about a few %. 

Besides the above-mentioned pure QED corrections another class of electro- 
weak higher order corrections should be included in the calculations. The 
gxcuw-f package has two schemes to do this: The running coupling constant 



and the scheme [p!6| . In the former the coupling constant of the fermion- 
fermion-Z vertex, Qf/z, is determined by evolving it from zero momentum 
transfer to the mass squared of the vv system, gf , which differs from one 
event to another. It varies according to the renormalization group equation 
(RGE) as 
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9ffz(q 2 z ) = g ff z{0) (i-£E Ctf log ^ 



-1 



1 



where a = 1/137.036 is the QED coupling at the zero momentum transfer, 
Ci the color factor, the electric charge, m 8 the mass of the z's fermion, 
and i runs over all massive fermions. In this scheme 9ffz(<lz) should be fixed 
event-by-event, but is common for all Z-exchange diagrams in the same event. 

The latter scheme is such that the weak couplings are determined through 
the weak-mixing angle, sin$w, which is given by 

• 2 Q rra{q 2 ) 1 

sin W 



v^M&l-Ar' 

where Mjy being the H^-boson mass and the muon decay constant. In this 
scheme the input is the precisely measured value of G^ at zero momentum 
transfer and a(q 2 ) is the QED coupling evolved to a typical energy scale of the 
processes which should be chosen by the users. In the present version of the 
program the loop corrections expressed by Ar in the formula is ignored(Ar = 
0). Then the difference between the running gffz scheme and the G^ scheme 
is to use the fixed energy scale(for G^) or the event-dependent energy scale(for 
running g ffZ ). 

The package grew-/ has a switch to choose either of the above corrections 
or no higher-order corrections. When users select the last option, the weak 
couplings are simply fixed by Mw and Mz through the on-shell relation, 
sin 2 ^ = 1 — ^pf , where M z is the mass of the Z-boson. 
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3 Numerical Results 



3.1 I7 test 

The numerical results for the single-photon emission are given in this section. 
All the three generators have been run and compared. The parameters used 
in the calculations are summarized in Table 1. 



M z 


91.187 GeV 


T z 


2.49 GeV 


M w 


80.22 GeV 




0.511 x 10- 3 GeV 


a(0) 


1/137.04 




1/128.07 



Table 1: Parameters used in the calculations. 

The weak coupling constant is determined scheme by scheme according 
to the selection of the higher-order corrections as explained in the end of the 
previous section. The visible photon is defined by Eq.(|l|) with E = 1 GeV 
and fl = 10°. 

First the results for are discussed. In this case four program packages 
are available: the analytic 0(a) calculations,^ NUNUGPV, the present one and 
K0RALZ. The last one also uses the complete matrix-elements for this process 
including the YFS exponentiation. The cross sections from four independent 
generators are visualized in Fig.3, and the numbers at four energy points 
are shown in Table 2. In these calculations a = 1/137.04 and sm 2 9w = 
1 — Mjy/M§ are commonly used for all four packages. The statistical errors 
to the numbers in the table are around 0.2% due to the numerical integration. 
The results from K0RALZ are very close to those from the 0(oi) calculation. 
Since the former includes the soft-photon exponentiation but the latter does 
not, one sees that the resultant difference in the total cross sections is less 
than 1%. In the grcz/z/7 package, however, QEDPS allows also the hard-photon 
exponentiation in the LL approximation. From the table this brings an effect 
of around 1.4%. By comparing three Monte Carlo program packages one can 
conclude that the theoretical uncertainty from the ISR corrections is about 
1%. The results of NUNUGPV are lower than grcz/z/7 by around 2%. 

2 The formulae given in Ref.H are used in this paper. 
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160 GeV 


170 GcV 


180 GeV 


190 GeV 


grcz/z/7 


2.470 


2.033 


1.709 


1.460 


O(a) corr. 


2.436 


2.006 


1.689 


1.449 


KORALZ 


2.437 


2.009 


1.697 


1.459 


NUNUGPV 


2.410 


1.987 


1.668 


1.430 



Table 2: Total cross sections of single-photon emission for the process in pb. The 
on-shell scheme without any higher-order electro-weak corrections is used for all packages. 

The differential cross sections with respect to the CMS energy of the pair 
isv, its longitudinal component, energy and transverse energy of the hard 
photon are compared in Fig.4. The shape of the distributions from these 
three packages (grcz/z/7, 0(a), and KORALZ) is in good agreement. A more 
detailed comparison between grcz/z/7 an d KORALZ on the same distributions 
is given in Fig.5, which shows the ratio of the results. Besides the overall 
factors these two are consistent. 

The effects of the higher-order corrections due to the running ffZ cou- 
pling are shown in Fig. 6 and Table 3 at four energy-points. The higher order 





160 GeV 


170 GeV 


180 GeV 


190 GeV 


grcz/z/7 only) 


2.470 


2.033 


1.709 


1.469 


t + running g ffz 


2.846 


2.334 


1.968 


1.678 


| for all neutrino species 


11.16 


9.70 


8.62 


7.85 



Table 3: Total cross sections of single-photon emission for the and all neutrino pro- 
cesses using grcz^7 in pb. 

effects from the running g^ z contribute to the total cross section by around 
13%. After summing up all the neutrino species the total cross sections 
are also summarized in the figure and the table. The contribution from the 
VF-exchange diagrams(last three diagrams in Fig.l) amounts about 25%. 

The scheme dependence of the higher-order corrections is found by grcz/z/7 
to be around 1% as seen from Fig. 7 and Table 4. Here q 2 = M"^ is chosen as 
the energy scale for the scheme. The discrepancy between two schemes 
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160 GeV 


170 GeV 


180 GeV 


190 GeV 


grcuwy + running gffz 


11.16 


9.70 


8.62 


7.85 


grcz/z/7 + G^ scheme 


10.98 


9.56 


8.53 


7.81 


NUNUGPV 


10.76 


9.41 


8.37 


7.69 



Table 4: Total cross sections of the single-photon emission for the all neutrino processes 
using grc^z/7 and NUNUGPV in pb. The latter package implements the G M scheme for the 
higher-order corrections. 

can be understood as follows. Two different energy scales are involved in this 
process: The colliding energy of the initial e + e~ and the Z-boson mass.f] The 
contributions to the cross section from these two energy regions are compati- 
ble in magnitude as seen from Fig.4. An expected difference in the total cross 
section could be estimated as 1 — a 2 (M%)/a 2 (El EP2 ) « 2% between these 
two scales. Hence the scheme which depends on a single fixed energy 
scale is not suitable for those processes that involve different energy scales. 

The differential distributions in G^ scheme obtained by grcz/z/7 are in 
good agreement with those by NUNUGPV as shown in Fig.8. 

3.2 27 test 

For the double-photon emission a similar comparison is also made between 
grcz/z/7 an d NUNUGPV. Both packages equip the exact matrix-elements for the 
double-photon emission associated with all neutrinos including z/ e z/ e (Fig.2) 
and the G^ scheme for the higher order corrections but different ways for 
the LL exponentiation for ISR, QEDPS or ^-depending structure function. 
Hence a difference occurs again in the treatment of the second configuration 
of the photon emission discussed in section 2. In applying grcz/z/7 to the 
double-photon emission, one has to take into account the possibility that one 
of the two photons supplied by the matrix-elements can be invisible (in Q ± ) 
and one photon from QEDPS becomes visible. The configuration in which both 
photons from the matrix-elements disappear into Q± is ignored safely. The 

3 There is another energy scale q 2 = for the real photon emission. This can be easily 
separated by using a fixed value of a — 1/137.04 at the e-e-7 vertex in the matrix elements 
and the ISR part. 
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total cross sections of the double-photon emission are summarized in Fig. 9 
and Table 5, where a slight difference is seen between two packages. The 





160 GeV 


170 GeV 


180 GeV 


190 GeV 


grcz/z/7 


0.687 


0.599 


0.537 


0.481 


NUNUGPV 


0.669 


0.578 


0.507 


0.472 



Table 5: Total cross sections of the double-photon emission for all neutrino processes 
using grcwy and NUNUGPV in pb. The latter package implements the G M scheme for the 
higher-order corrections. 

differential distributions are in good agreement with each other as shown in 
Fig.10. 

4 Anomalous WW^f Coupling 

In the grcuwy package the anomalous coupling of the W-W-7 vertex is pre- 
pared. The program includes only those terms which conserve C and P 
invariance, derived from the following effective Lagrangian||17||: 

L eff = -ze[(l + A gij ){WlW» - W^W^A" + (1 + Ak^W^A^ 

~ A/f2 ^A* ^ J' 
m W 

where W^ v = d p W u - d u W^ A^ u = d^A v - d u A^ Here Ag ly , Ak 7 and A 7 
stand for the anomalous coupling parameters which vanish in the standard 
model. 

To find the region sensitive to the anomalous TGC in the photon phase- 
space, the energy and angular distributions are examined at Ag l7 = 0, A/t 7 = 
— 10 and A 7 = 0. This is shown in Fig. 11. The sensitive region, as expected 
from a glance at the diagrams, locates in the high energy side and large angle 
region, except for around the 'Radiative Return', where the annihilation 
dominates the process. Then the experimental cuts of 20 GeV< _E 7 <55 
GeV or 65 GeV< _E 7 and 45° < # 7 < 135° are effective to enhance the signal 
from the anomalous TGC. The total cross section as a function of Ak 7 and A 7 
with the experimental cuts at the CMS energy at 170 GeV is shown in Fig. 12. 



10 



Once 200pb~ 1 of the accumulated luminosity and 100% of the acceptance are 
assumed, 128 ± 11 events will be observed. The error given here is only a 
statistical one. The sensitivity of —9 to +6 for A 7 and —3 to +10 for Ak 7 
is expected at the three-sigma limit if one supposes only the statistical error 
of the measurements. Since the higher order corrections also change the 
total cross sections under the above cuts, as pointed out in Ref . || , a careful 
investigation is needed on the systematics from those corrections. 

5 Conclusions 

An event generator grew-/ for the processes e + e~ — > uuj and e + e~ — > 
uujj has been proposed. The exact matrix elements for single- and double- 
photon emission are prepared by GRACE. QEDPS is used as a tool for the ISR 
corrections. For the case the total cross sections and the hard-photon 
distributions of grevvy are compared with those from the 0(a) calculation, 
KORALZ and NUNUGPV. It is found that the theoretical uncertainty for the 
ISR corrections is under control at the 1% level. The systematics of the 
scheme from the double energy scales involved in the reaction is estimated to 
be around 1%. The energy spectrum of the hard-photons is in a reasonable 
agreement with independent calculations from KORALZ and NUNUGPV up to 
the double-photon emission. 

Concerning v e a similar comparison with NUNUGPV has been done, though 
in lacking 0(a) calculations. With this process there is some opportunity 
to measure the anomalous TGC at the W-W-j vertex which appears in the 
single energetic events at the LEP-2 when the luminosity is enough accumu- 
lated. Since only the experimental observable is the total cross section, a 
careful investigation of the systematic errors must be done. 
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Figure 1: Feynman diagrams of the process e+e — > u e u e j. Only the first two diagrams 
appear in and v T neutrino production. 
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Figure 2: Feynman diagrams of the process e+e — > v e v e ^. Only the first six diagrams 
appear in ty, and v T neutrino production. 
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Figure 3: Total cross sections of the single-photon emission for the process in pb. 
The on-shell scheme without any higher-order electro-weak correction is used for all the 
program packages. 
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Figure 4: Differential cross sections with respect to the CMS energy of vv system, its 
longitudinal component, energy and transverse energy of the hard photon obtained by 
grevvj (solid histograms), 0(a) calculations (stars), and K0RALZ (circles) at the CMS 
energy of 170 GeV. 
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Figure 5: Ratio between and KORALZ for the same distributions as in Fig. 4. The 

error bars show the statistical errors of the numerical integration. 
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Figure 6: Total cross sections of the single-photon emission for the and all neutrino 
processes using grcvvj in pb. 
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Figure 7: Total cross sections of the single-photon emission for all the neutrino processes 
using grcvvy and NUNUGPV in pb. The latter implements the scheme for the higher 
order corrections. 
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Figure 8: Differential cross sections with respect to the CMS energy of the vv system, 
its longitudinal component, energy and transverse energy of the hard photon obtained by 
grci^7 (solid histograms) and NUNUGPV (circles) at the CMS energy of 170 GeV. 
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Figure 9: Total cross sections of the double-photon emission for all the neutrino processes 
using grci^7 and NUNUGPV in pb. The latter implements the G M scheme for the higher 
order corrections. 
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Figure 10: Differential cross sections with respect to the CMS energy of the vv system, 
its longitudinal component, energy and transverse energy of the hard photon obtained by 
grcvv^ (solid histograms) and NUNUGPV (circles) at the CMS energy of 170 GeV. 
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Figure 11: Energy and angular distributions of the hard photon obtained by grcvvj 
with the standard TGC (solid histograms) and the anomalous TGC with Ak 7 = — 10 at 
the CMS energy of 170 GeV. The distributions of the ratio between the standard and the 
anomalous TGC are also shown. 
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Fi gure 12: Total cross sections of the single hard-photon events as a function of A 7 
(left) and Ak 7 (right) at the CMS energy of 170 GeV. The dashed lines are the expected 
experimental upper limit at three standard-deviations of the statistical error from the 
standard model parameters after 200pb~ 1 data accumulation. 
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Appendix 



A.l How to install the program 

The source code is available by anonymous ftp from ftp.kek.jp in the 
directory kek/minami/grcnna. The grcz/z/7 system contains the complete set 
of FORTRAN sources for e + e~ — ■> z/z/7 and e + e~ — > 2/2/77 matrix-elements 
and the four libraries, i.e., BASES/SPRING, CHANEL, QEDPS and utilities for 
kinematics. Those source codes are written in FORTRAN77. grcz/z/7 nas 
been developed on HP-UX, but should run on any UNIX platform with a 
FORTRAN compiler. 

The procedure of installation is as follows: 

1. Editing Makefile. 

The following macros in Makefile should be taken care of by users 
themselves. For example, GRCNNADIR defines the directory name where 
grcz/z/7 is installed. The values of FC and FOPT define the relevant 
compiler name and option for your system. The other macros can be 
left as they are. 

GRCNNADIR = directory where grcz/z/7 are installed. 
LIBDIR = directory where libraries are installed. 

(default is $ (GRCNNADIR) /lib.) 
BINDIR = directory where an executable is installed. 

(default is $ (GRCNNADIR) /bin.) 
MACHINE = [hpux|hiux|sgi|dec|sun] 
FC = FORTRAN compiler command name. 

FOPT = FORTRAN compiler options. 

CERNLIBS = CERNLIB including the jetset library. 

For HP-UX, -L/cern/pro/lib -1 jetset 74 

-lpacklib -lkernlib -L/lib/pal . 1/ -lm 

2. Compilation. 

By executing command make install the executable of the interface 
program is generated at BINDIR. Furthermore four libraries, BASES/SPRING, 
CHANEL, QEDPS and kinematics utility library, are generated in LIBDIR. 
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A. 2 How to run the program 

According to the number of hard-photons, two directories are prepared: 
$ (GRCNNADIR) / la and $ (GRCNNADIR) /2a. Three FORTRAN files and a 
Makefile are stored in each directory. 

For the numerical integration, BASES: 

• mainbs.f: A main program of the numerical integration. 

• usrprm.f: A subroutine for parameter setting. 
For the event generation, SPRING: 

• mainsp.f : A main program of the event generation. Event-loop is also 
included. 

The subroutine, usrprm.f, is used commonly for both of BASES and SPRING. 

Users can set the parameters of the calculations by editing the file usrprm . f . 
The key words are as follows; 

* CM energy 

EO = 170. dO 

* Hard Photon condition 

Ea = 1 . dO 
THEa = 10. dO 

* channel selection 

* Ineut=l : sum up all neutrino species 

* =2 ; nu-e only 

* =3 : nu-mu (nu-tau )only 
Ineut = 1 

* QED correction ********************************************* 

* QEDPS switch ! 

* lps=0 : No ISR ! 

* =1 : ISR by QEDPS ! 
Ips=l ! 

* QED correction ********************************************* 
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* Higher order correction scheme ***************************** 

* Running alpha ! 

* lruna=0 : On-shell scheme ; fixed alpha ! 

* =1 : On-shell scheme ; running alpha ! 

* =3 : Gmu scheme ! 
Iruna=l ! 

* I 

* Running width ! 

* lrunw=0 : fixed width ! 

* =1 : running width ! 
Irunw=l ! 

* Set QED alpha ! 

alphaO= 1.0d0/137.0359895d0 ! 

if (Iruna.ne . 1) then ! 

alpha = 1.0d0/128.07d0 ! 

else ! 

alpha = 1.0d0/137.0359895d0 ! 

end if ! 

* Higher order correction scheme ***************************** 

* WWgamma coupling 
* 

* angla : gl 

* andka : delta-k 

* anlma : lambda 
angla =1 
andka =0 
anlma =0 

The parameters for the anomalous TGC are available only for the single- 
photon case. 

Users can proceed with the calculations as follows: 
i) Change directory by typing 

°/„ cd $ (GRCNNADIR) / la 



=1 (SM) 
=0 (SM) 
=0 (SM) 
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or 

°/„ cd $ (GRCNNADIR) /2a 

ii) Edit usrprm . f to set user parameters. 

iii) Create an executable integ for the integration by typing 

°/„ make integ 

iv) Numerical integration is actually performed by typing 
°/„ integ 

The results of the integration step are displayed on the console as well 
as written in an output file bases .result. The total cross section 
in pb and the estimated statistical error are shown on the last line, 
under Cumulative Result, in the table of the Convergence Behavior 
for the Integration step. The differential cross sections are also 
printed as a function of the energy, scattering angle of each particle and 
invariant masses of any two final particles. The probability distribution 
of the integrand is written in a file bases. data which will be used in 
the event generation step by spring. 

v) Before running the event generation, users may edit mainsp.f to set 
additional parameters if needed and call the user's own analysis rou- 
tines. 

The following is the structure of the generated mainsp.f, where four- 
momentum of all particles including soft-photons generated by QESPS 
are stored in the common/lujets/ in the JETSET format when subpro- 
gram sp21nd is called in the event-loop: 

Program mainsp 
implicit real*8(a-h, o-z) 
external func 



real*4 p,v 
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common/lujets/n,k(4000,5) ,p(4000,5) ,v(4000,5) 



mxtry = 50 
mxevnt = 1000 
do 100 nevnt = 1 , mxevnt 

call spring ( func, mxtry ) 



* 

call sp21nd 

* 

* ============================================= 

* ( user_analysis based on the common lujets ) 

* ============================================= 

* 

100 continue 



stop 
end 

iv) Create an executable spring for event generation by typing 

°/„ make spring 
vii) Start the event generation by typing 

% spring 

Information concerning the event generation will be written in the 
spring. result file. Users should pay attention to the histograms 
generated in this step. The distributions of the generated events are 
superimposed with the character "0" on the histograms generated in 
the integration step. These two distributions should be consistent with 
each other within the statistical error of the generation. For the details 
of the output files of BASES and SPRING, users can consult Ref.[[L5|. 
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